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O Abstract 

We derive general and explicit expressions for the unrenormalized and renormalized dressed 
propagators of fermions in parity-nonconserving theories with inter-generation mixing. The mass 
^■i eigenvalues, the corresponding mass counterterms, and the effect of inter-generation mixing on 



their determination are discussed. Invoking the Aoki-Hioki-Kawabe-Konuma-Muta renormaliza- 
tion conditions and employing a number of very useful relations from Matrix Algebra, we show 
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explicitly that the renormalized dressed propagators satisfy important physical properties. 
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I. INTRODUCTION 



The aim of this paper is to derive general and explicit expressions for the unrenormalized 
and renormalized dressed propagators of fermions in parity-nonconserving theories with 
inter-generation mixing, and to discuss their important physical properties and implications. 

The results presented here immediately apply to the Standard Theory of Elementary 
Particle Physics, usually referred to as the Standard Model (SM), as well as its extensions. 
As has been known for a long time, the quark fields are subject to inter-generation mixing, 
as implemented by the Cabibbo-Kobayashi-Maskawa (CKM) [1] quark mixing matrix. Since 
neutrino oscillations have been observed experimentally and lower mass bounds have been 
established, the lepton fields are known to also undergo inter-generation mixing. An early 
treatment of flavor- changing self-energies, both for leptons and quarks in bound states, 
which, however, focuses on finite renormalization effects, may be found in Ref. [2|. On the 
other hand, our treatment is quite general and takes into account the full mixing amplitudes. 
The renormalization of the CKM matrix has been recently discussed by several authors; see, 
for example, Ref. |3[ and references cited therein. Mixing renormalization has also been 
worked out for theories involving Majorana neutrinos 4|. 

This paper is organized as follows. Section [III discusses the derivation of the unrenormal- 
ized dressed propagators. The mass eigenvalues, the corresponding mass counterterms, and 
the effect of inter-generation mixing on their determination are also analyzed. Section IIHI 
discusses the renormalization of the dressed propagators. Invoking the Aoki-Hioki-Kawabe- 
Konuma-Muta (AHKKM) renormalization conditions and employing very useful relations 
from Matrix Algebra, it is shown explicitly that the renormalized dressed propagators satisfy 
important physical properties. Section IIVI contains our conclusions. The Appendix explains 
how to derive the two- loop expression for the mass eigenvalues presented in Sec. [Til an d how 
to express the mass counterterms in terms of the unrenormalized self-energies. 



II. UNRENORMALIZED DRESSED PROPAGATOR OF MIXED FERMION 
SYSTEM 

As is well known, the unrenormalized mass matrix can be brought to diagonal form 
with non-negative eigenvalues by means of bi-unitary transformations on the left- and right- 
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handed fields. On this basis, the unrenormalized inverse propagator is where 

Iii{#) = ty-rr$)5 ij -V ij {t), (1) 

i, j are flavor indices 1 and the self-energies are given by 

E y (0 = [f{B + ) l3 + [A + ) i3 ] a + + [f{B^ + (A_) y ] a_. (2) 

In Eq. (j2J), (B±)ij are Lorentz-invariant functions of p 2 and a± = (1 ± 7s)/2 are the 

chiral projectors. 2 

Equations ([1]) and can be written in compact form, as 



-T + )a + + (^_-T_)a_, (3) 

where S + and T + are matrices defined by 

(S±)ij = - (B ± ) iv (T±)^ = m°<% + (At)^. (4) 

The unrenormalized dressed propagator is iP{jf) = i(/(^)) -1 . 3 Writing (/(^i))' 1 = (fU+ — 
V + )a + + (fU- — V_)a_, we find the relations 

S + V + + T_U + = 0, (5) 

+ T + E7_ = 0, (6) 

p 2 S + U^+T^ = 1, (7) 

p 2 S^U+ + T+V+ = 1, (8) 

where 1 stands for the unit matrix. 

In order to express U± and V± in terms of S± and T±, we first solve for V_ in Eq. ([H]) and 
insert the result in Eq. (j7J). This leads to 

EL = [j?S + -T-{SJ)- 1 T+]- 1 , (9) 

K_ = -(5_)- 1 T + E/_. (10) 



1 In this paper, repeated indices are not summed, unless a summation symbol is explicitly included. 

2 Throughout this paper, we adopt the notational conventions of Bjorken and Drell |5|. 

3 Here and in the following, the matrix iP{f) is referred to as the unrenormalized propagator. The particle 
propagators are the elements of this matrix, namely £Py(ii). An analogous denomination is used in Sec. lIIII 
for the renormalized propagator iP{p) and the renormalized particle propagators iPijiii). 
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Next we solve for V + in Eq. (JSJ) and insert the result in Eq. (jSj), which leads to 

U + = [p 2 5_-T + (5 + )- 1 T_]- 1 , (11) 
V+ = -{S + )~ 1 T.U + . (12) 

More convenient forms for U± are obtained by writing 

U- = [{p 2 -T-iS-r'T+iS+r 1 ) S+]- 1 = (5 + )- 1 (p 2 - DC)-\ (13) 
U + = [(p 2 -T + (S + )- 1 T4S_)- 1 )SJ\- 1 = (S„)- 1 (p 2 -CD)-\ (14) 

where 

C = T+(S+)- 1 , D = T_(5_)- 1 . (15) 
It is also convenient to introduce the matrices 

E = (S+y 1 ^, F = (S-.)- 1 ^. (16) 

Using Eqs fl9|)-f lT4|) . f lT5|) . and f|T6|) . the unrenormalized dressed propagator is given by iP, 
where 

P = {f + £?)( i S_)- 1 (p 2 - CL>)- 1 a + + (f + F)(S + )- 1 (p 2 - DC)- x a_, (17) 

which is fully expressed in terms of the self-energy matrices S± and T±. The matrices 
(p 2 — CD)^ 1 and (p — DC)~ l are related by similarity transformations, as 

(p 2 - CD)- 1 = C(p 2 - DC)- 1 ^ 1 = D-\p 2 - DC)~ X D. (18) 

Writing 

^- GD ^ = ^-CD Y V-^' W-DC) ' (19) 
where a + and a_ are the corresponding adjoint matrices, 4 we see that the determinants are 

equal and that a + and a_ are related by the same similarity transformations as in Eq. (|T8l) . 

Thus, the squared mass eigenvalues M 2 are the zeros of det(p 2 — CD), namely they satisfy 

det(M 2 - Y(M 2 )) = 0, (20) 
Y{p 2 ) = (CD)(p 2 ). (21) 



4 Given a square matrix M, in this paper the adjoint matrix Adj M means the transpose of the matrix 
whose elements are the cofactors of M (see, for example, Ref. Q.) We recall that the cofactor Cy of the 
element rriij of M is (— iy +: > times the determinant of the matrix obtained by deleting the i-th row and 
the j-th column of M. 
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The off-diagonal elements of Y(p 2 ) arise from inter-generation mixing and are, therefore, of 
0(g 2 ) or higher, where g is a generic weak-interaction gauge coupling. As a consequence, 
if terms of 0(g A ) are neglected, only the diagonal elements of p 2 — Y(p 2 ) contribute to the 
determinant, and the eigenvalues are of the form 

M 2 = Y rt (M 2 ) + (D(g*), (22) 

where Yu(p 2 ) denotes Yu(p 2 ) in the absence of 0(g 4 ) contributions. If, instead, terms of 
0(g 4 ) are retained, but three-loop contributions and higher are neglected, there are two ad- 
ditional effects: (i) there are now terms of 0(g 4 ) in Yu(p 2 ) and (ii) the non-diagonal elements 
Yij(p 2 ) (i 7^ j) contribute to the determinant. As a consequence, the mass eigenvalues are 
now of the form 

M} = Y U (M 2 ) + £ + 0{g\ g"a s ). (23) 

In the Appendix, we outline the derivation of Eq. ([231 and show how Eqs. (jl]), (fT5l) . (l2Tj) . and 
( f2~3l can be used to express the mass counterterms in terms of the unrenormalized self-energy 
functions A ± and B ± , in the approximation of neglecting three- loop-contributions. 



III. RENORMALIZED DRESSED PROPAGATOR OF MIXED FERMION SYS- 
TEM 

In order to renormalize P [cf. Eq. (fTTj) ]. we recall that the unrenormalized propagator 
is the Fourier transform of (0| T(\l/ (x)\l/ (0)) |0), where the zero superscripts denote the 
unrenormalized fields. In our case, they are column and row fields with components labeled 
by flavor indices. In the following discussion, we assume for simplicity that the fermions are 
stable. Decomposing the fields into right- and left-handed components, as 

= + ^ W° = ^+ + ^° , (24) 

where = a±^ and = ^a T , and taking into account the effect of the chiral pro- 
jectors a±, it is easy to see that the first, second, third, and fourth terms of P arise from 
(0| T(ty°y°_) |0), (0| T(*°.#°.) |0), (0| T(^° |0), and (0| T(^°J° + ) |0), respectively. 
Shifting the fields according to 

^ + = Z^ + , *°=* + (zf)t, (25) 

= z y^_ 5 y°_ = y_(z 1 ! 2 )\ (26) 
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where \l/± are the renormalized fields, we see that the four terms in P are multiplied on 
the left and right by various combinations of z]± and (Z± 2 y factors. Since the time- 
ordered products are now expressed in terms of renormalized fields, in order to obtain the 
renormalized propagator iP, we must divide out such factors. Specifically, the first term in 
P must be multiplied on the left by Z_^ 2 and on the right by (Z^ 2 y , the second term by 
Z^ 2 on the left and (Z_^ 2 y on the right, the third term by Z + ^ 2 on the left and (Z + ^ 2 y 
on the right, and the fourth term by Z_^ 2 on the left and )t on the right. 

Thus, the renormalized propagator is iP, where 

p = (pZZ 1/2 + Z- 1/2 E)(S-)-\p 2 -CD)- 1 (ZZ 1/ y a+ 

+ (fZ- 1/2 + ZZ 1/2 F)(S + )-\p 2 - DC)- x (Z; 1/2 )ta_. (27) 

Recalling Eqs. (Tl5|) and (fl6|) . we see that the third and fourth terms are related to the first 
and second ones, respectively, by the exchange + -H- — . 

We now note that the Z' 1 / 2 factors in Eq. ([27)1 can be absorbed in a redefinition of the 
self-energy matrices S± and T±, namely 

S± = {Zfys+Zf, T± = {Zfyr ± zf. (28) 

Using Eq. ( |28l) . P can be written in the compact form 

P = {f + S)(^_)-V - CD)- 1 a+ + {f + F)(5+)-V - t)d)- x a., (29) 

where 

c = f + sz\ D = f_SZ 1 , E = 5; 1 f_, F = SZ 1 f + . (30) 
In particular, CD and CD are related by a similarity transformation, as 

CD = {zHyCD{zty- 1 , (31) 

so that det(p 2 — CD) = det(p 2 — CD) and the mass eigenvalues are the zeros of either 
determinant. The matrices S±, T±, C, D, E, and F are the renormalized counterparts of 
S±, T±, C, D, E, and F, respectively. 

In analogy with Eq. (TlBl . we have the relations 

( P 2 - cd)- 1 = c(p 2 - dc)-^- 1 = b-\p 2 - bcy 1 b. (32) 
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We note that CD and FE are also related by a similarity transformation and so are (p 2 — 
CD)- 1 and (p 2 - FE)- 1 , namely 

S^CDS- = FE, Sz\p 2 - CD)- 1 ^ = {p 2 - FE)- 1 . (33) 
Interchanging + <H- — , we obtain 

S- 1 DCS + = EF, S- 1 {p 2 -DC)- 1 S + = {p 2 -EF)-\ (34) 

Using Eqs. (ED}, ([32}, (EH, and §S>, Eq. ([29]) can be cast in the alternative form 

P = a_(p 2 - FE)- l SZ l {f + C) + a + (p 2 - EF^ST 1 ^ + D). (35) 

It differs from Eq. ( 129]) in that the chiral projectors a± are on the left side of the expression. 
In both Eqs. (f29j) and (|35|) . the cofactors of a_ and a + are related by the exchange + •<-)■—. 
Writing 

(p 2 - CD)- 1 = (p 2 - DC)' 1 = ^— (36) 

det(p 2 -OD) det(p 2 -DC) 

(p 2 - FE)' 1 = (p 2 - EF)- 1 = ^— — (37) 

det{p 2 -FE) det{p 2 -EF) 

where a± and j3± are the corresponding adjoint matrices (cf. Footnote H]), the similarity 
relations in Eqs. (E2])-(E4]) tell us that 

det(p 2 - CD) = det(p 2 - DC) = det(p 2 - FE) = det(p 2 - EF), (38) 

and 

S-'a ± = P ± S-\ (39) 

a relation that plays an important role in our discussion of the propagator's properties. We 
recall that in the previous equations, S±, T±, C, D, E, F, a±, and f3± are functions of p 2 . 

We now turn our attention to the renormalization conditions. As emphasized in the semi- 
nal work of AHKKM 7;j, a fundamental physical property of the renormalized propagator iP 
is that, as f — > m n , where m n is one of the mass eigenvalues, the pole (fl — m„) _1 should be 
present only in the diagonal element iP nn of the renormalized propagator matrix. In order 
to implement this property, as well as the conventional requirement that the pole residue 
equals the imaginary unit, AHKKM proposed suitable conditions on the renormalized in- 
verse propagators, which were described both graphically and mathematically. 



Recalling Eq. fl3]), in our general matrix notation the renormalized inverse propagator is 
—il(]f), where 

I{f) = (A - T+)a+ + {fS_ - f_)a_. (40) 
An alternative expression is 

I{f) = a_{fS + - f_) + a + (pS„ - f+), (41) 

where the chiral projectors a± are placed on the left. The homogeneous AHKKM renormal- 
ization conditions read 

u n (jf>)Ini{f) = 0, (42) 
iln{f)Un{f) = 0, (43) 

where u n {jj>) is a spinor that satisfies fu n {i>) = m n u n (fi), u n {j/>) its hermitian adjoint, and n 
and / are flavor indices. 

Inserting Eq. (1401) into Eq. ()42l) . we have 

mj ± (m 2 n )-f±(m 2 n )} =0. (44) 

J nl 

Multiplying on the right by S± (Trvfyij, summing over /, and remembering the definitions in 
Eq. (130|) . this becomes 

C n j(m 2 n ) = D nj (m 2 n ) = m n 5 nj , (45) 

which implies 

(CD) nn (ml) = ml, (CD) nj (m 2 n ) = (j^n), (46) 

with the analogous result for (Z)(7)(m^). 

Inserting Eq. (14T1) into Eq. (1431) . recalling the definitions in Eq. ( |30l) . and carrying out 
the analogous analysis, we obtain 

E in (m 2 n ) = F in (m 2 n ) = m n S in , (47) 

which leads to 

(EF) nn (m 2 n ) = ml (EF) m (m 2 n ) = (i ^ n), (48) 

and the analogous result for (F£')(m^). 

Equation (146]) tells us that, as p 2 — > m 2 n , all the elements in the n-th row of p 2 — CD and 
p 2 — DC vanish. Therefore, the only non- vanishing cofactors of (p 2 — CD) and (p 2 — DC) are 
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those corresponding to the elements of that row, namely the cofactors C n \. Since the adjoint 
matrices are the transpose of the cofactor matrices (cf. Footnote H]) , we conclude that the 
only non- vanishing elements of a±(m 2 l ) are those in the n-th column, namely the elements 
(a±) in (m 2 l ). Similarly, from Eq. (|48p we see that, as p 2 — > m 2 , all the elements in the n-th 
column of p 2 — FE and p 2 — EF vanish. Consequently, the only non- vanishing elements of 
$±{ m n) are those in the n-th row, namely {I3±) n j{'m 2 l ). In combination with Eq. ( )39|) . these 
results imply that, as p 2 — > m 2 , the only non-vanishing elements of the matrices S^ 1 a± and 
P±S^ are the diagonal nn elements (S^ 1 a±) nn (m J 2 ^ = (fi±SZ l ) nn (rn%) . Thus, 

(S-'&iUml) = 0±S-%(m 2 n ) = (z or j ? n). (49) 

To examine the effect of these results on the renormalized propagators, we insert Eq. f )36|) 
and Eq. (l37|l into Eq. ( 129]) and Eq. (!35|) . respectively. Recalling Eq. (|38|) . we obtain 

p = if + E)(S-)- 1 &+a + + {ft + F)(g + )- 1 d_a_ 
detQo 2 - CD) 

from Eq. (|29|) . and 

p = aJ + {S-)-\f + C) + aJ-{S + )-\f + D) 
det(p 2 -FE) 

from Eq. (1351) . 

Using Eqs. (|4"5]) . (|47|) . and ( 149|) . one readily verifies that, as p 2 —> m 2 , the only non- 
vanishing elements in the numerators of Eqs. fj50|) and ( l51~j) are, in fact, the diagonal nn ele- 
ments. Thus, the explicit expressions of the renormalized propagator iP, given in Eqs. (|29|) , 
( |35|) . ( 1501) . and ( loTl) . indeed satisfy the fundamental physical property that the (ft — m n ) _1 
pole is present only in the diagonal element iP nn of the propagator matrix. 

The inhomogeneous AHKKM renormalization conditions are 

-Inn{ft)Un{ft) = U n {ft) , (52) 



— m 



n 



Un(ft)Ln(ft)— = U n (ft). (53) 

ft Win 

Inserting Eq. ( l4"Uj) into Eq. (I5"2"j) . expanding the numerator about ft = m n , and using Eq. 
we find the renormalization conditions 



5_(m n ) + m n m n (S + + 5_) - T + - T_ 

(S-WK) = (5 + )„ n (m 2 ), (51) 



where the prime symbol stands for the derivative with respect to p 2 , evaluated at p 2 
Inserting Eq. (jHJ) into Eq. fl53l) . we obtain the same result. 

In order to analyze the effect of Eq. (I5"3j) . we evaluate the residue of the (fl — m n ) _1 pole 
in P using Eq. floDl) and focus on the a + term. We expand det(p 2 — CD) about p 2 = w? n 
through 0(p 2 — m 2 )). Since p 2 = w? n is a zero of the determinant, the first term vanishes, 
and we have 

det(p 2 - CD) = \det(p 2 - CD)] ' {p 2 - m 2 n ) + • ■ ■ . (55) 
Using the well known expression 

(det M)' = Tr (M' Adj M) , (56) 

the r.h.s. of Eq. (155]) becomes Tr j&^m 2 ) 1 — (CD)' | (p 2 — m 2 ) + ■ ■ • . Multiplying by 
f — m n , taking the limit f — > m n , and recalling Eqs. ( )47l) and ()49l) . we see that the residue 
of the (p — rrin)^ 1 pole in the a + term of Eq. ( 1501) is 

(5'_ 1 Q; + ) rtn 



Res_i 



Tr < a 



1 - (CD)' 



(57) 



Here and in the following, it is understood that all the functions are evaluated at p 2 = vr? n . 
To simplify this expression, we insert S-Sl 1 = 1 in the argument of the trace. Recalling 
again Eq. ( H9l) . we find 



Tr 

and the residue becomes 



1 - (CD)' 



Res i 



S'_ 1 q;_ 



1 - (CD)' 



5_ 



{[l - (CD)'] S\ 

L L J J nn 

Taking into account Eqs. ( H5|) - p71) . Eq. (|59|) becomes 

1 



(58) 



(59) 



Resj 



m n (S + + S-) -T+-T- 



(60) 



S- + m. 

Thus, the renormalization condition of Eq. ( 15^1) indeed implies that 

Resj_ = 1. 



(61) 



Calling Res_ the residue of the 
analysis shows that 



m n ) 1 pole in the a_ term of Eq. (l5Dj) . an analogous 
Res_ = 1. (62) 
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We conclude that, when the inhomogeneous renormalization condition of Eq. (1521) is imposed, 
the poles in our explicit expressions for the renormalized propagator [cf. Eqs. ff29|) . f )35|) . fl50|) . 
and (l5"Tj) ] have residues i. 

IV. CONCLUSIONS 

We derived general and explicit expressions for the unrenormalized and renormalized 
dressed propagators of fermions in parity-nonconserving theories with inter-generation mix- 
ing [cf. Eqs. f|T7|) . f[2"§j) . (I3l)j) . floU]) . and (1ST])]. We analyzed the determination of the mass 
eigenvalues and the corresponding mass counterterms in the approximation of neglecting 
three-loop contributions [cf. Eqs. (|23|) and ( 1A9|) ]. In particular, we discussed the effect of 
inter-generation mixing on these determinations. Using the AHKKM renormalization con- 
ditions and applying very useful relations from Matrix Algebra, we showed explicitly that 
our renormalized dressed propagator [cf. Eqs. ( 1291) . ( |35i) . ( |50|) . and (15T1)]. which is valid to all 
orders in perturbation theory, satisfies important physical properties. In turn, this demon- 
strates in a clear manner that the AHKKM renormalization conditions are also valid to any 
order of perturbation theory. 
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Appendix A: Appendix 

In this appendix, we outline the derivation of Eq. ( 12 3 p in the approximation of neglecting 
three-loop contributions and show how it can be applied to express the mass counterterms 
in terms of the basic self-energy functions {A±)ij and {B±)ij in Eq. (T5]). For simplicity, we 
consider the three-generation case. 

As explained in the paragraph containing Eqs. (120]) and (T2T]) . the mass eigenvalues are 
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the zeros of det(p 2 — Y (p 2 )), where Y(p 2 ) = (CD)(p 2 ) and the matrices C and D are defined 
in Eq. (PJ. Using Eqs. ©, (EES]), and we find 

y = (M ) 2 + Z, (Al) 

where M° is the diagonal bare mass matrix with elements and 

Z = (M°) 2 £-(l + BJ) + M°(A_ + A_B_ + + M°B + (1 + 5+)!/° 

+ + B+)M° + A + M°B_ + M°B+M°B_ + (A2) 

In Eq. (lAip . we have separated out the squared bare mass term (M ) 2 and the one- and 
two-loop contributions contained in Z. We recall that A±, B±, Y, and, consequently, Z are 
functions of p 2 . It is further convenient to split 

(M ) 2 = M 2 + 5M 2 , (A3) 

where M is the renormalized mass matrix whose elements are the mass eigenvalues and 5M 2 
is the mass counterterm matrix. Thus, 

Y = M 2 + X, (A4) 

where 

X = 5M 2 + Z. (A5) 

We note that, in Eq. (IA4I) . M 2 contains the zeroth-order terms, while X contains the one- 
and two-loop contributions. 

Neglecting three-loop contributions, in the three-generation case the eigenvalue equation 
det(p 2 — Y(p 2 )) = becomes 

(p 2 -Yn)(p 2 -Y 22 )(p 2 -Y 33 )-(p 2 -Y n )Y 23 Y 32 -(p 2 -Y 22 )Y 13 Y 31 -(p 2 -Y 33 )Y 12 Y 21 = 0. (A6) 

Consider the neighborhood of p 2 = M 2 , where Mi is one of the mass eigenvalues: dividing 
by {p 2 - Y 22 )(p 2 - Y 33 ), we have 

*23*32 



(AT) 

P — 1 33 P r 22 



(P 2 -^22)(P 2 -V33) 

The factors I23I32, ^13^31, and Y\ 2 Y 2 \ are of two-loop order or higher. As p 2 — > M 2 , we see 
from Eq. (lA4j) that, to leading order, we have p 2 — y 22 = — Mf and p 2 — Y33 = M 2 — M|. 
Thus, neglecting three- loop contributions, as p 2 — » M 2 , Eq. (1A7I) reduces to 

M 2_ Y (m 1 o^xm?) (y 13 y 3 i)(M 2 ) 
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which is a particular case of Eq. (|23|) . 

Recalling Eqs. fl23|) . (lAlj) . and flA3j) . the mass counterterms are then 



Mf - ill 2 



^)-E ( ^ )( 5 } > (A9) 



Mf-M- 



where Z is defined in Eq. ( 1A2|) . In the last equality of Eq. ( 1A9I) . we have replaced Yy — )■ Z^, 
since both are equal when % ^ j [cf. Eq. C IA 1 j) ] ■ 

We note that, subject to our approximation, the amplitudes involving linear powers of 
A± and B± in Eq. (1A2|) contain both one- and two-loop contributions. 

Using Eq. f ]A2j) . we find for the diagonal terms 



Z u = {m° i f{B + + B_ + B 2 + + B 2 _) u + m°(A + + A_ + A + B + + A_B_ + B + AJ) U + (A + AJ) U 

3 

+ J2 H(A + UB-h + rr%m!j(B + ) ij {B-.) ji ] . (A10) 
i=i 

We note that depends not only on the bare fermion masses m" and m°- displayed in 
Eq. flAlOp . but also on additional ones present in the loop diagrams. We refer generically 
to the latter as mf. Consistently with our approximation, in the contributions of two- loop 
order, we replace the bare masses m", m", and by the mass eigenvalues M iy Mj, and Mi, 
respectively. In the contributions of one-loop order, we replace 



Ml-Zf{Ml)X'\ (All) 



and similarly for m°. In Eq. flAllj) . the superscript (1) stands for the one-loop contribution, 
namely 

Zf = M 2 (B^ + B { l ] )(Mf) + Mi{A^ + AL 1} )(Mf), (A12) 

with an analogous expression for Z^\ 

The contributions involving ZijZji with j ' ^ i in Eq. (|A9[) are already of two-loop order or 
higher, so that in the off-diagonal amplitudes Z^ with j ^ i, we simply replace m-, m", — > 
Mi, Mj, Mi. In this way, subject to the approximation of neglecting three-loop contributions, 
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the mass counterterms 5Mf given in Eq. (1A9I) are fully expressed in terms of the basic self- 
energies A±(M?) and B±(M?) of Eq. (j2J) and the mass eigenvalues. 
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